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Recent scaling results for the AC conductivity of ionic 
glasses by Roling et al. [Phys. Rev. Lett. 78, 2160 (1997)] 
and Sidebottom [Phys. Rev. Lett. 82, 3653 (1999)] are 
discussed. It is shown that Sidebottom's version of scaling is 
completely general. A new analytical approximation to the 
universal AC conductivity of hopping in the extreme disorder 
limit, the "diffusion cluster approximation," is presented and 
compared to simulations and experiments. 



Disordered solids have AC electrical properties remark- 
ably in common. These solids, in fact, have so similar 
frequency-dependent conductivity a(uj) that ionic con- 
duction cannot be distinguished from electronic. Even 
the temperature-dependence of <t(w) is "quasi-universal." 

The class of disordered solids with quasi-universal AC 
behavior is large, including polycrystalline and amor- 
phous semiconductors, organic semiconductors, ionic 
conductive glasses, ionic melts, non-stoichiometric crys- 
tals, ionic or electronically conducting polymers, metal 
cluster compounds, transition metal oxides, and macro- 
scopic mixtures of differently conducting phases like 
organic- inorganic composites, etc. Each class contains 
hundreds of different solids and there is a huge literature 
on their AC conductivities. 

It is possible for almost all disordered solids to scale AC 
data at different temperatures into one single curve. This 
so-called master curve gives the dimensionless AC con- 
ductivity a = cr(cj)/a(0) as function of a scaled dimen- 
sionless frequency a). The existence of a master curve is 
sometimes referred to as the "time-temperature superpo- 
sition principle" (TTSP). The master curves of different 
solids - while not identical - are surprisingly similar. This 
quasi- universality was recognized gradually in the 1970's 
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The common AC features of disordered solids are 
the following 0]: At low frequencies conductivity is 
frequency-independent. Around the dielectric loss peak 
frequency to m [|| AC conduction sets in, and for uj 3> ui m 
cr(u)) is close to a frequency power-law. The exponent 
is between 0.7 and 1.0 [Q. As temperature is lowered 
the exponent goes to 1.0. In a log- log plot the AC con- 
ductivity is much less temperature-dependent than the 
DC conductivity. A final ubiquitous observation is the 
Barton-Nakajima-Namikawa (BNN) relation |^-[Tl|] con- 
necting dielectric loss strength Ae |6| , dielectric loss peak 
frequency, and DC conductivity - er(0) = p Ae eo 0J m - 
where p is a numerical constant close to one. 



To construct the master-curve frequency must be 
scaled by u> m . Because the dielectric loss strength is only 
weakly temperature-dependent while cr(0) and uo m are 
both Arrhenius, the BNN-relation implies oj m ~ cr(0). 
Thus, the existence of a master curve is conveniently 
summarized into 
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where C may depend on variables like charge carrier con- 
centration n, temperature T, high frequency dielectric 
constant, etc. 

Recently, there has been renewed interest in scaling 
and universality of AC data for ionic glasses 12-13] (ex- 



cellent reviews of glass ionic conduction have been given 
by Nowick and coworkers Jll| and Ngai JIt]]). In 1997 
Roling, Happe, Funke, and Ingram showed that C oc n/T 
for sodium borate glasses Jl2|. This year, however, Side- 
bottom showed that in general scaling is not achieved by 
C oc n/T - instead three quite different ionic conductive 
systems all obey the following scaling relation jl5| , 
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The purpose of this paper is three-fold. First, we put 
the recent scaling results into a historical perspective. 
Second, we show that Eq. (0) is generally valid. Finally, 
Sidebottom's experimental master curves are compared 
to two models, the symmetric hopping model and the 
macroscopic model. In this connection a new analytical 
approximation to the universal AC hopping conductivity 
is presented. 

First, the historical perspective. In a series of papers 
towards the end of the 1950's Taylor analyzed the di- 
electric properties of ionic glasses in accordance with the 
Debye equation with a spread of relaxation times p8[ . 
He showed that the dielectric loss M for all glasses fell 
on a single plot against scaled frequency. In 1962 Isard 
relabeled Taylor's axis by plotting dielectric loss against 
log of the product of frequency and resistivity and thus 
essentially arrived at AC scaling in the form given in 
Eq- (0) [^9)- Since then Eq. (|l|), which we shall refer to 
as "Taylor-Isard scaling," has been used in several dif- 
ferent contexts. For instance, Taylor-Isard scaling was 
used by Scher and Lax in 1973 in their famous papers 
introducing the continuous time random walk approxi- 
mation pO[ j, by Summerfield and coworkers in 1985 for 
amorphous semiconductors ^J,^2|, by van Staveren and 
coworkers in 1991 for metal-cluster compounds pq ], and 
by Kahnt the same year for ionic glasses (24] . 
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We now proceed to comment on the scaling princi- 
ple of Roling and coworkers, C oc n/T in Eq. ([!]), and 
Sidebottom's recent interpretation of its occasional vi- 
olations. If q is charge, fjj jump rate, and d jump 
length, Sidebottom bases his arguments on the expres- 
sions ct(0) oc nq 2 d 2 f H /k B T and Ae oc nq 2 d 2 /k B T ||. 
Clearly, to obtain Cj frequency should be scaled by fn, 
thus leading to C oc n/T if all other parameters are con- 
stant. But since it is reasonable to assume the jump 
length would increase as concentration decreases, Side- 
bottom argues that one cannot expect C oc n; on the 
other hand the n-dependence of d is taken care of by Eq. 
(^) in which the unknown jump length is eliminated by 
scaling with the measured Ae. 

In our opinion Sidebottom's arguments are largely cor- 
rect, but Eq. (|^) is much more general than it appears 
from his Letter. Consider hopping of completely non- 
interacting charge carriers on a cubic lattice. In this 
model conductivity is n times charge carrier mobility 
and for given lattice jump frequencies the mobility scales 
with d 2 . This shows that C oc nd 2 as Sidebottom has 
it. However, if the lattice is homogeneous with just one 
jump frequency conductivity is frequency- independent; 
in order to have strongly frequency-dependent conduc- 
tivity lattice jump frequencies must vary many decades. 
In such more realistic cases the expressions used by Side- 
bottom for cr(0) and Ae do not apply - there simply is 
no unique ju- Still, we find in our simulations that these 
more realistic hopping models always obey Eq. (||). This 
leads to the question: When is Eq. (^) obeyed? 

Equation (Q) apparently expresses too informations: 1) 
TTSP is obeyed, and 2) the scaled frequency is given by 
Cj = [eoAe/cr(0)] u>. We now show that 1) mathematically 
implies 2): TTSP implies the existence of some function 
cr(Cj) where Cj is the scaled frequency. Expanding to first 
order in Cj leads to a = 1 + iACj (where A is real because 
c*(u>) = a{—uj)). Since a = aa(0) we have a = a(0) + 
iAoj(j(0). On the other hand, from the definition of Ae 
H one has a = er(0) + iuiAeeg for cj — > 0. Equating the 
two expressions for a leads to Q = A^ 1 [Aeeo/cr(0)] oj, 
i.e., Eq. (p|) |25[| (the numerical value of A is irrelevant). 
Thus, Eq. (§)follows from TTSP alone. 

We next compare the data discussed by Sidebottom to 
two models, the macroscopic model and the symmetric 
hopping model. The macroscopic model considers AC 
conduction in a random mixture of phases with different 
(frequency-independent) conductivities |2^,^]. The AC 
conductivity is obtained by solving Maxwell's equations. 
Computer simulations have shown p^ , p^] that in the ex- 
treme disorder limit (where the local conductivities vary 
over many decades) the AC conductivity is accurately 
described by the effective medium approximation (EMA) 
universality equation |2(||27],^9| , 

cr In a = iuj . (3) 

The symmetric hopping model [p3-p2[ , which models AC 



conduction via random walks of non-interacting parti- 
cles on a cubic lattice with random symmetric Arrhenius 
nearest- neighbor jump rates (joe~^ E , where E is the ran- 
dom energy barrier) , also exhibits universality in the ex- 
treme disorder limit. Even for this model does the EMA 
lead to Eq. (||) Q. For hopping, however, Eq. (||) does 
not give an accurate representation of the universal &{Cj) 
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FIG. 1. Numerical results for the symmetric hopping 
model in 3 dimensions with periodic boundary conditions |m| . 
Reported results are averages over 100 independent N x N x N 
cubic lattices, a) Real part of a(ui) for the Box distribution 
of energy barriers (p(E) = 1, < E < 1), scaled according 
to Eq. (§); (3 = 20 (N = 14), f3 = 40 (N = 24), f3 = 80 
(TV = 32), and (3 = 160 (N = 64). As increases the data 
converge to a universal curve, b) The apparent exponent 
/i = d In a' /d In u plotted versus a . Data are shown for 5 dif- 
ferent energy barrier distributions, p(E), see |^7| for details. 
/3 = 160 (N = 64hwhere @ = P/p(E c ) and E c is the "perco- 
lation energy" |j27|]33| . The universal curve is independent of 
the energy barrier distribution. Estimated values of Aeeo/P 
for different distributions agree within 1%. The numerical 
data are compared to three analytical approximations: EMA 
(Eq. (§)), PPA ( §5)), and DCA (Eq. §). The universal curve 
lies between EMA and PPA and is well approximated by DCA 
with do = 1.35. 
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Before comparing model predictions to experimental 
data we present results from new simulations of the sym- 
metric hopping model. The simulations were done on 3- 
dimensional samples with periodic boundary conditions, 
using a method based on exactly solving the master equa- 
tion p4| . Figure la shows a 1 - the real part of a - for the 
Box distribution of energy barriers. The frequency axis 
is scaled according to Eq. (||). Clearly, as (3 = 1/fcgT in- 
creases a' converges to a single curve. As seen in Fig. lb 
showing the apparent exponent /x (a' ~ uj^) as function 
of a' ', the universal curve is the same for different energy 
barrier distributions. Figure lb compares the simulations 
to three analytical approximations, the EMA universal- 
ity equation Eq. , the percolation path approximation 
(PPA) [^,^5), and a new "diffusion cluster approxima- 
tion" (DC A), which we describe now: 

The reason for AC universality in hopping is the fact 
that in the extreme disorder limit conduction takes place 
on the "percolation cluster" formed by the links with 
largest jump rates until percolation ^,^,^|; the PPA- 
idea (27j]3j| is to regard the conducting paths on the per- 
colation cluster as strictly one-dimensional. The phys- 
ical idea behind EMA is to replace the inhomogeneous 
lattice by an effective homogeneous medium determined 
self-consistently |3(],|3l]]. Our hopping simulations show 
that the truth is somewhere between PPA and EMA 
p4j , somewhere between one-dimensionality and homo- 
geneity. Most likely, this is because both approxima- 
tions ignore the fact that conduction takes place on some 
complex subset of the percolation cluster. This "diffu- 
sion cluster" must be smaller than the backbone (de- 
fined by removing dead-ends of the percolation cluster, 
fractal dimension=1.7 |54 38 1) and larger than the set 



of red bonds (those that, when cut, stops the current, 
fractal dimension=l.l J34|,|38|]). At present this is all we 
know about the diffusion cluster and its dimension do 
is regarded below as a fitting parameter. To derive the 
DCA equation we use EMA in do dimensions. In the 
extreme disorder limit EMA implies (3^]3j| Iiict oc sG 
where, if p(k) = j- J2i=i cos(fej), sG is the following 
integral sG = f ^ u . ,^\ d - — „ , %u} u — t^tt- Whenever 

1 < d < 2 one finds sG oc (iuj/a) do ^ at relevant Q 
frequencies. Thus, after rescaling frequency we arrive at 
the DCA equation, 



Inff = I - 
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As is clear from Fig. lb, the solution of this equation 
for do = 1.35 gives an excellent fit to the universal AC 
hopping conductivity. 
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FIG. 2. Experimental data for AC conduction of 3 ionic 
systems scaled according to Eq. (0) (f = lj/2tt): 

a) (Na 2 0) x (Ge0 2 )i- x (Sidebottom ||) 

b) 0.4Ca(NO 3 ) 2 - O.6KNO3 (Howell et. al. 

c) (K 2 S) x (B 2 S 3 )i-x (Patel 

Experimental data are compared to the macroscopic model 
(represented by EMA), the universal AC conductivity of the 
symmetric hopping model in the extreme disorder limit (nu- 
merical data), and DCA (Eq. (jij)). (DCA is scaled to agree 
with the numerical data for the hopping model, since it cannot 
be scaled by using Eq. (|) (Ae = 00).) 
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In Fig. 2 the data discussed by Sidebottom jTa] 
are compared to the macroscopic model represented by 
the EMA universality equation (|J) (dashed line), to 
the hopping model simulations in the extreme disor- 
der limit (filled squares), and to the DCA (full line, 
Eq. (0)). Figure 2a shows Sidebottom's own data for 
(Na20) x (Ge02)i-x glasses, Fig. 2b shows data for the 
ionic melt 0.4Ca(NO 3 ) 2 - O.6KNO3 measured by Howell 
and coworkers p^ , pO| , and Fig. 2c shows Patel's data 
for the thioborate glass system (K 2 S) x (B 2 S3)i_ x Jl5| , p| . 
In all three systems we find a frequency range where ex- 
perimental data agree well with the symmetric hopping 
model. 

The symmetric hopping model is an extremely sim- 
ple model which contains no fitting parameters but still 
gives a good fit to experimental data as shown above. 
The model does not include Coulomb interactions and 
it allows an arbitrary number of charge carriers at each 
site. We note, however, that the model is mathematically 
equivalent to that obtained by linearizing (with respect to 
the electric field) a hopping model with energy disorder 
and Fermi-statistics |$(J. The macroscopic model does 
not work very well for the ionic systems in Fig. 2. This 
is interesting because the macroscopic model via Gauss' 
law includes all effects of Coulomb interactions. Appar- 
ently, the continuum description behind this model does 
not appropriately reflect the actual microscopic disorder. 

To summarize, we have briefly outlined the history of 
AC scaling in disordered solids and showed that the scal- 
ing version Eq. (^) must be obeyed whenever conductiv- 
ity scaling is possible at all (TTSP obeyed). We have 
presented a new analytical approximation to the univer- 
sal AC conductivity of hopping in the extreme disorder 
limit, the diffusion cluster approximation (DCA), and 
shown that it gives an excellent fit to simulations 
Also, DCA (and thereby the extreme disorder limit of 
the symmetric hopping model) agrees well with the ionic 
data discussed by Sidebottom [|l5| , in a frequency range 
depending on the system. 
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